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ABSTRACT 

Mid-infrared (IR) observations of polycyclic aromatic hydrocarbons (PAHs) and molecular hydrogen emis- 
sion are a potentially powerful tool to derive physical properties of dense environments irradiated by intense 
UV fields. We present new, spatially resolved, Spitzer mid-IR spectroscopy of the high UV-field and dense 
photodissocation region (PDR) around Monoceros R2, the closest ultracompact H 11 region, revealing the spa- 
tial structure of ionized gas, PAHs and H2 emissions. Using a PDR model and PAH emission feature fitting 
algorithm, we build a comprehensive picture of the physical conditions prevailing in the region. We show 
that the combination of the measurement of PAH ionization fraction and of the ratio between the H2 0-0 S(3) 
and S(2) line intensities, respectively at 9.7 and 12.3 yum, allows to derive the fundamental parameters driv- 
ing the PDR: temperature, density and UV radiation field when they fall in the ranges T = 250 - 1500 K, 
Hh = 10"^ - lO^cm"^, Go = 10^ - 10^ respectively. These mid-IR spectral tracers thus provide a tool to probe 
the similar but unresolved UV-illuminated surface of protoplanetary disks or the nuclei of starburst galaxies. 
Subject headings: ISM: lines and bands — ISM: molecules — infrared: ISM 



1. INTRODUCTION 

Dense (> 10"^ cm"^) and high UV field (> 10"^ times the 
standard interstellar radiation field in units of the Habing field, 
written Go hereafter) photodissociation regions rule the en- 
ergy balance, and thus evolution, of some of the most fun- 
damental astrophysical objects such as protoplanetary disks, 
planetary nebulae and starburst galaxies. One of the best 
tools to probe this UV illuminated matter is spectroscopy 
in the mid-infrared (5-15 yum, hereafter mid-IR), because it 
provides information on both the gas and small dust grain 
properties (polycyclic aromatic hydrocarbons and very small 
grains, hereafter PAHs and VSGs) . Ideally, one would like 
to spatially resolve the emission of these components in order 
to study their variations as the UV field is attenuated. Un- 
fortunately, such observations are very hard to achieve (for 
the moment) on one hand because large ground-based tele- 
scopes, providing arcsecond angular resolution, are restrained 
in wavelength coverage, while on the other hand, space borne 
telescopes have diameters that are usually too small to resolve 
the sources. Because it is the closest ultra compact H 11 re- 
gion at a distance of 850 p c, Monoceros R2 (Mon R2, see 
IWood & Churchwell I [T989I and iHoward etaH Il994h consti- 
tutes one of the rare exceptions where one can resolve the 
PDR between the Hii region and molecular cloud. In this 
Letter, we present and analyze the mid-IR PAH and molecu- 
lar hydrogen emissions in Mon R2 based on spatially resolved 
Spitzer spectral mapping. 

2. OBSERVATIONS 

Mon R2 was observed using the Infrared Spectrograph 
(IRS) onboard Spitzer, in the low resolution mode (A/ A A = 
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60 - 127) as part of the "SPECHII" program (PI C. Joblin). 
The data were obtained in the spectral mapping mode. The 
full spectral cubes of the SLl and SL2 modules were built us- 
ing the CUBISM software ( Smith et al.ll200 7) from the basic 
calibrated files retrieved from the Spitzer archive (version S19 
pipeline). The two cubes were then assembled to provide the 
full SL cube of 26x36 positions in space and -170 points in 
wavelength, ranging from 5 to 14.5 yum. 

3. OBSERVATIONAL RESULTS 

An overview of our observations is presented in Figdl The 
ionized gas in the H 11 region, traced by the [Ne 11] line at 12.8 
yum, appears confined in a spherical region, as seen b y radio 
continuum observations (IWood & Church welTI [198 9^ and in 
higher angul ar and spectral r esoluti on [Neii] g round based 
observations (iTakahashi et al.1 120001: iJaffe et al.1 12003.) . The 
cometary shape of the H 11 region is well seen in our [Ne 11] 
map and peaks at about 0.12 erg s"^ cm"^ sr"^ in the sur- 
roundings of the B 1 star IRS 1 . The fact that [Ne 11] maxium 
contours do not include the exact position of IRS 1 is due to 
the saturation of the IRS detectors at this position, implying 
that the intensity coul d not be measure d, though likely peak- 
ing there as shown by ITakahashi et al.l (2000) and Jaff'e e t al.l 
(2003). PAH emission is present everywhere in the region, 
but for clarity in Fig[T] we only present the region where it 
is the brightest in the 11.3 jim band (1 to 4 xlO"^ erg s"^ 
cm"^sr"^) forming a filamentary/shell structure that surrounds 
the Hii region. The H2 0-0 S(3) rotational line at 9.7 jim 
peaks in a filament which lies between the Hii region and 
the cold and dense molecular gas traced by CS J=5-4 (Fig.[T]) 
transition. The intensity of the H2 0-0 S(3) line is of the or- 
der of 1-4 xlO-4 erg s"^ cm'^ sr'^ The H2 0-0 S(2) line 
at 12.3 ju m follows a similar spatial distribution. As in the 
Orion bar (iTielens et al. 1993), the spatial distribution of these 
H2 lines is not correlated with the PAH emission contrary 
to what is seen in lower UV field PDRs like the Horsehead 
nebula (Habart et al. 2003; Compiesne et al. 2007) or the p- 
Ophiucus filament ( Habart et al 2005). In the following we 
investigate the origin of these structures from a physical point 
of view, using PAH emission and modeling the H2 excita- 
tion in the PDR. To simplify the task we select three diff'erent 
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Fig. 1. — Overview of the Mon R2 region seen with Spitzer. Left: Observed IRS spectra for PDR 1, 2 and 3 (continuous Hne) and fit by the model (diamonds, 
see Sect l4.2.l1 . Middle: The blue contours represent the intensity of the [Neii] line (0.02 to 0.1 erg s~^ cm~^sr~^ in linear steps), the red contours the intensity of 
the PAH 1 1.3 jum band (1.8 to 3.4 xlO~^ erg s~^ cm~^sr~^ in linear steps) and the green contours the intensity of the H2 0-0 emission in the S(3) rotational line at 
9.7 /um (1.5 to 4.5 xlO~^ erg s~^ cm~^sr~^ in linear steps). The background map shows the CS /=5-4 emission presented in Choi et al. ( 2000). Low emission is 
in black. Numbers indicate the positions of the selected PDRs. Right: Cut in the map along the yellow box, showing the stratified emission of the diff'erent lines. 



zones named PDRl, PDR2, and PDRS as displayed in Fig[T] 
PDRl was chosen to be representative of the region where the 
maximum of H2 emission is found (PAH emission is weaker 
but present). PDR2 is the transition from the H 11 region to the 
molecular gas. At these two positions it is well seen that H2 
gas is found further from the star than the PAH emission (see 
cut in FigdJ. Finally, we positioned PDRS in a region situated 
further from IRS 1 which is a smoother transition from ionized 
to neutral and molecular gas. For comparison, we also con- 
sider the mid-IR spectrum obtained integrating the whole IRS 
cube over the area imaged. 

4. PROBING THE PHYSICAL CONDITIONS OF A HIGH UV FIELD, 
HIGH DENSITY PDR 

4.1. H2 emission as a probe of gas density and temperature 

In dense and high-UV PDRs like Mon R2, the excita- 
tion of the H2 p ure rotational lines is dominated by in- 
elastic collisions (Le Bourlot et al.|[T999h . To quantitatively 
probe the role of density nu and radiation field Go on the 
excitation of the lowest energy rotatio nal levels of H2, we 
use the revised Meudon PDR code (iLe Petit et al.l 120061: 
iGoicoechea & Le Bourlot Il2007l) . for a large grid of radiation 
fields and densities above 10^ cm"^. Since the critica l den- 
sity of the S(S) line is ~ 5 x 10^ cm'^ (iLe Bourlot etanil999l) 
the lowest energy rotational levels are thermalized when nn > 
5x10"^ cm"^. As a consequence the S(S)/S(2) line ratio scales 
with the gas temperature and Trot - T, where Trot is the 
rotational temperature of the associated rotational transition. 
In high UV high density PDRs, the main mechanism heat- 
ing the gas is the photo-electric heating by electrons ejected 
from very small dust grains and PAHs. Thus, the photoelectric 
heating efficiency will depend on the ability to eject electrons 
from the grain surface (less efficient as grains become posi- 
tively charged) and on the density of electrons in the gas with 
which charged grains recombine and neutralize. In PDRs, 
low energy electrons are provided by the ionization of carbon 
atoms, and as essentially all the carbon is ionized, the electron 
density depends on gas density and carbon abundance. Over- 
all, as the elemental abundace [C]/[H] is constant, an increase 
of the gas density (thus of electron density) reduces the grains 
charge and increases the photoelectric heating efficiency (i.e., 
the gas temperature). This explains the dependence of the 
S(S)/S(2) ratio with density seen in Fig 121 

4.2. PAHs as a probe of radiation field 
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Fig. 2. — Results from the Meudon PDR code for high UV field/density 
conditions, a) Variations of the H2 0-0 S(3) absolute intensity (colors and 
contours) as a function of the hydrogen nuclei density nn and intensity of 
the radiation field Gq. b) Variations of the S(3)/S(2) line ratio (colors and 
contours) as a function of nn and Gq. Circles indicate the intersection be- 
tween observed values of 75(3) (panel a)) and lso)lh{2) (panel b)) and the 
Go derived using PAH ionization fraction for PDR 1,2 and 3. 



4.2.1. Mid-IR emission model 

In order to analyse t he full Spitzer mid -IR c ube we use 
a model adapted from "Jobli n et aH (l2008h and IB erne et aH 
(|20Q9t) . This emission model includes the four major dust 
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Fig. 3. — Contours of PAH*^ and PAH+ intensity (resp. green and blue) 
overlay ed on [Ne ii] intensity map, built from the results of the fit obtained 
for each spectrum of the cube of Mon R2. Positions of IRS 1,2,3 marked as 
in Fig.[T] 



components, introduced in these previous works: VSGs 
(small carbonaceous dust grains), free neutral PAHs, posi- 
tively ionized PAHs (PAH^), and finally PAH^ that are large 
(with Nc > 100 carbon atoms) charged PAHs. We decom- 
pose the observed mid-IR em ission using the 4 template spec- 
tra presented in Berne et al.' ("2009). In addition, we con- 
sider continuum emission that is simply modelled with two 
slopes. Furthermore, we include the eff'ect of extinction 
by multiplying the mid-IR spectrum by (1 - ^"^^)/t^ where 
= Cext(^) • Ni, CextW being the extinction c ross s ection 
per nucleon, taken from 'Weig artner & Drainel (1200 1|) with 
Ry = 5.5 and Ni, the column density on material on the line 
of sight is left as a free parameter in the fit. The values Ay 
of visual extinction corresponding to the derived Ni (related 
by Ni = 1.8 X lO^^Ay) for each PDR, and used to correct 
the measured H2 lines intensities as a function of their wave- 
length, are given in Table [T] We fit all the spectra of the cube 
using this technique, thus providing the spatial distributions 
of the emission of each component. An example of these fits 
is provided in Fig. [T] and spatial distributions of PAH^ and 
PAH^ populations are presented in Fig. [3] 

4.2.2. Estimating Gq using PAH ionization ratio 

As shown by mo dels (!Tielens jr2005i) and ob servations (e.g. 
Bregman & Temi 1 12005 : Flage v et al.l 120061: iGalliano et al] 
2008D, the PAH ionization ratio ([PAH^]/[PAH% depends 
on the parameter Go ^fT /rie where T is the gas temperature 



and He the electron density. Ionization of PAHs will influ- 
ence the emission cross section of these molecules, so that 
the ratio between the 6.2 to 11.3 jim bands ih.ilhi.?) will 
increase with increasing Go ^ /rie. Assuming that most elec- 
trons are provided by the ionization of carbon, we can write 
thai He - x(C^)nH - [C]/[H]nH, with [C]/[H] th e elemental 
abund ance of carbon assumed to be 1.6 x 10 ""^ (ISofia et al. | 
I2OO 4I) . Finally, using the empirical law of Galliano et al.l 
©08), we can relate the h.i/hi.d to Go(T /lO^^^ /hh by: 



Go(T/loY^/nH - (1990[C]/[if])x((/6.2//ii.3)-0.26) (1) 

Using the ratio between the 6.2 and 1 1.3 /urn bands for neutra l 
and ionized PAHs found in Table 1 in Rapacioli et al.l (|2005|) 
we relate h.ilhi.?, to the ionized to neutral PAH density ratio 
^^1^, using Eq. (2) in lJoblin et al.l ([T996h : 



/6.2//11.3 = 1.12 X 



1 + 1^x0.85 
l + I^xO.29' 



(2) 



5. APPLICATION TO MON R2 PDRS 



The high value of 1^(3) we observe require both Go > 
1 X 10^ and uh > 1 X 10^ cm"^ (Fig. 2 and ' Burton etall 
1 19921: [Kaufman et al.h2006) . This implies that the observed 
Is (3) /Is (2) ratio, after correction for extinction (Table [B, al- 
lows to derive the approximate densities of PDR 1,2,3 in 
Fig|2l-b. In addition, since these lines are thermalized, the 
gas temperature, T, is coincident with the rotational temper- 
ature. Trot, and we can infer T for the diff'erent PDRs (Table 
[21). To calculate the rotational temperature we have assumed 
an ortho-to-para ratio of 3, that is the equilibrium value for 
temperatures larger than 100 K. Finally, using the value of 
Go(T / 10^)^ ^^/hh estimated with the PAH ionization fraction 
(or h.2lhi3), and the above derived T and hh we can de- 
rive the intensity of the radiation field that illuminates the 
three PDRs and position precisely PDRs 1,2 and 3 in Fig. 
2. The found values for hh, T, Go, (Table ^ are consistent 
with Mon R2 being a dense and highly U V irradiated PDRs 
as estimated from othe r molecular lines (IChoi et al.[ [20001: 
[Rizzo et aH [2003L [2005*). The found values for density are 
lower than in these previous works, because the present trac- 
ers (PAHs and H2) probe the outermost cloud regions directly 
exposed to the UV radiation field that are usually warmer 
(>500 K) and less dense than regions situated deeper into the 
molecular cloud (see cut in Fig. 1). Finally, the estimated pa- 
rameters using this technique, for the entire Mon R2 region 
(see Table 2), are more consistent with the ones derived for 
spatially resolved PDRs 2 and 3. 

6. CONCLUSIONS 

The main results we have presented in this Letter are: (/) 
Mon R2 is a unique template of high-UV/high-density PDRs 
that is spatially resolved, similar to the Orion bar (where the 
UV field is slightly lower), (//) our observations are consis- 
tent with PDR models which predict that in high-UV and 
dense PDRs, the excitation of pure rotational H2 lines is col- 
lisional and depends on gas temperature, (///) for this rea- 
son, the spatial distributions of PAHs and warm H2 are dif- 
ferent, contrary to what is seen in cool PDRs, (iv) however, 
because PAH emission is present in the extended region illu- 
minated by the radiation, their ionization fraction can be used 
to probe the intensity of the UV radiation field, even where 
H2 emission peaks. Thus, the mid-IR spectrum of dense and 
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TABLE 1 

Observational diagnostics in PDRl, PDR2 and PDR3. Observational errors are given for H2 lines. 



H2 lines (corrected from extinction) fit of PAH bands 





^S{2) 




ho)lh{2) 


Av 


h.ilhi.z 


[PAH+] 


GQ^|T|\0^ 




(erg cm-^ s-^ sr'h 




(K) 








(k1/2 cm~3) 


PDRl 


2.3±0.15 X iQ-4 


9.6±o i X 10-4 


4_2±o.30 


16 


1.1 


0.03 


0.28 


PDR2 


1.6±o i X 10-4 


2.1±o i X 10-4 


13+0.15 


12 


1.8 


1.60 


0.49 


PDR3 


1.8±o.i xlO-4 


2.0^0.1 X 10-4 


1 1+0.15 


19 


2.0 


2.44 


0.56 


Mon R2^ 


8.0±0-i X 10-5 


10.3±0-i X 10-5 


1 3±o.io 


9 


1.9 


1.9 


0.52 



spectrum obtained integrating the IRS cube over area imaged 



TABLE 2 

Derived physical parameters for PDRl, PDR2 and PDR3. Errors are the 

RESULT OF PROPAGATION FROM THOSE IN TaBLE 1 . 





UH (cm-^) 


TrotO^) 


Go 


PDRl 


4.3*0-5 X io5 


574*25 


1.6*0-2x105 


PDR2 


4.0*0-4 X io4 


331+19 
314 


3.3*0-3 X io4 


PDR3 


3.7*0-3 X io4 


3.7*0-2 X io4 


Mon R2 


4.0*0-4 X 104 


321+18 


3.7*0-4 X io4 



highly irradiated PDRs appears as an efficient probe of the 
physical conditions in these environments. The ratio between 
the molecular H2 0-0 S(3)/S(2) line intensities allows to di- 
rectly compute the density and temperature of the gas, while 
the measurement of the ionization fraction of PAHs allows to 
derive the intensity of the radiation field. We show that the 
derived parameters for the entire Mon R2 region are consis- 
tent with those found for the spatially resolved PDRs 2 and 3. 
This is likely because the dense PDRl occupies only a small 
part of the whole region, and therefore dilution eff'ects imply 
that the spatially averaged spectrum is dominated by emission 
from lower density PDRs. Thus, we suggest that this spectral 
methodology could be useful to derive the dominant physical 
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conditions in PDRs that are not spatially resolved in the mid- 
IR. In particular, the PDRs at the surface of proto planetary 
disks around Herbig Ae/B e stars (iBerne et al."2009') or in the 
inner rim of T-Tauri disks (lAgundez et al. 2008) are expected 
to be the targets where such an analysis can be applied. The 
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by the emission of PDRs having similar conditions as those 
described in this paper (see Fuent e"etan[200 8 for the case of 
M82). Thus, the present methodology could be a useful tool to 
derive global properties of galaxies, in connection with mas- 
sive star formation activity, using forthcoming James Webb 
and SPICA space telescopes that will observe these emission 
features at low and high redshifts. 
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